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fluctuation-dissipation regime 
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We show that the observed non-monotonic behavior of the thermal conductance between two 
nanoparticles when they are brought into contact is originated by an intricate phase space dynamics. 
Here it is assumed that this dynamics results from the thermally activated jumping through a rough 
energy landscape. A hierarchy of relaxation times plays the key role in the description of this 
complex phase space behaviour. Our theory enables us to analyze the heat transfer just before and 
at the moment of contact. 
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1. Introduction. — The knowledge as to how heat 
transfer takes place at interfaces between different mate- 
rials in nanoscale systems is an issue which has received 
considerable attention recently. This phenomenon, in- 
volved in high-precision lab techniques as well as in the 
field of active nanotechnology research P, [l], is poorly 
understood and the theoretical bases behind it are as of 
yet sketchy. Despite the difficulty to demonstrate the 
near-field enhancement of radiative heat transfer, many 
experiments have been performed in this area such as 
the archetypical of the heat transfer between a nano-tip 
and a surface [3] . As a consequence of the near-field ra- 
diative interactions, transitions among the energy states 
of the nanoparticles (NPs) related with structural rear- 
rangements occur. This phenomenon is similar to the 
Forster energy transfer 0]. 

The current literature on the subject of radiative en- 
ergy exchange at the nanoscale is implicitly or explic- 
itly based on the validity of the fluctuation-dissipation 
theorem Q involving fluctuating electric dipoles. In the 
particular case of two NPs which are brought together, 
the near-field radiative energy transfer between them de- 
pends greatly on their characteristic size D and separa- 
tion distance d. When the space between them is small, 
i.e. d around 2D, the energy exchange can be consid- 
ered as due to multipolar interactions in the fluctuation- 
dissipation regime [6] related to exponential or Debye 
relaxation. Nonetheless, when the NPs get even closer 
to each other up to contact d near to D, the inter- 
action energy becomes of the same order of magnitude 
as the unperturbed energy of the NPs which lose any 
crystalline structure they may have. This suggests that 
the NPs are trapped in a rough energy landscape sub- 
jected to a kind of glassy behavior common to a great 
variety of complex systems such as glasses 0] and pro- 
teins (H, Q, just to mention some examples. In such situ- 
ations a general methodology is lacking. Characterizing 
the energy landscape in terms of an order parameter, 
the heat exchange process is related to the diffusion of 



non-interacting quasi-particles through the order param- 
eter space. Besides, systems displaying glassy behavior 
have a multitude of metastable states that they visit after 
overcoming an irregular distribution of energy barriers. 
These systems exhibit aging and memory effects inher- 
ent to an activated dynamics characterized by a hierarchy 
of relaxation times TO]. Therefore, the NPs cannot be 
treated as a thermodynamic system at equilibrium and 
consequently, the fluctuation-dissipation theorem cannot 
be invoked to explain the energy transfer nor is a mul- 
tipolar expansion applicable. Then, the question arises 
as to how to understand the heat exchange between two 
bodies separated by a nanometric distance in such ex- 
treme and nevertheless common cases in experiments. 

In this Letter, we have provided a general theory to 
describe the radiative heat exchange in a general sit- 
uation. The observed non-monotonic behavior of the 
thermal conductance between two NPs when they are 
brought into contact, does suggests an activated process 
frequently found in a wide range of situations correspond- 
ing to complex systems [Tl|. The nontrivial topology of 
their phase space is inherent to the rough energy land- 
scape arising from the internanoparticle interaction. In 
the framework of mesoscopic nonequilibrium thermody- 
namics we flnd an expression of the radiative heat 
flux and derive the heat conductance which depends on 
the Fourier transform of a hierarchy of relaxation times. 
We then compare our result with previous molecular dy- 
namic data reproducing the strong enhancement and 
subsequent decay of the thermal conductance observed 
when the NPs are in contact. In our approach the heat 
exchange is interpreted as a current of quasi-particles 
seen as stationary waves bounded by both NPs which 
gives a result in good agreement with the simulation. 

2. Mesoscopic nonequilibrium thermodynamic analy- 
sis. — We will follow here the classical approach of the ac- 
tivated dynamics in which the system diffuses by thermal 
hopping through a rough energy landscape. The land- 
scape structures are determined by the strength of the in- 
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teractions between the microscopic constituents and with 
external agents, both global and local. 

Heat exchange by near-field radiation is based on 
Coulomb interaction. As a consequence of this interac- 
tion, transitions occur among the energy states of the 
NPs related with structural rearrangements. Hence, let 
us assume that the energy exchange process is given by a 
diffusion current of quanta of energy, quasi-particles [13| 
emitted and absorbed by bodies at different tempera- 
tures. The irreversible process represented by the exis- 
tence of this current can be analyzed in the framework 
of mesoscopic nonequilibrium thermodynamics based on 
the assumption of the validity of the second law in the 
phase space. Hence, let us consider a nonequilibrium gas 
of quanta characterized by the probability density p(r,i), 
where T = (p, x) and p, x are the momentum and posi- 
tion of a quasi-particle, resp ectively. According to the 
Gibbs entropy postulate 14 , [l^ , the density functional 
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is the nonequilibrium entropy of the system, where S^q. 
is the equilibrium entropy of the gas plus the thermal 
bath and p^q, (F) is the equilibrium probability density 
function. Changes in the entropy are related to changes 
in the probability density which, since the probability is 
conserved, are given through the continuity equation 
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where d/dV — {d/dp,\7), V = d/dx. The continuity 
equation ^ defines the probability current J {T,t) = 
{JxjJp) to be determined after computing the entropy 
production which follows by combining the time deriva- 
tive of Eq. ^ and © to give 
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This is an equation which according to the second law 
and analogous to the form of the entropy production in 
macroscopic diffusion [l^ , expresses the entropy produc- 
tion in terms of the nonnegative product of a thermody- 
namic current J (r,i) and its conjugated thermodynamic 
force d/dT [fj,{T,t)/T] , the gradient of the nonequilibrium 
chemical potential 
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which derives directly from Eq. ([T]). The physical mean- 
ing of the nonequilibrium chemical potential Q becomes 
clearer when we rewrite Eq. ([1]) 
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which is the Gibb's equation in the phase space. Since 
—T5S = 6A^ with A the nonequilibrium free energy, 
the nonequilibrium chemical potential (jl]) can be inter- 
preted as the free energy per unit of probability mass. 
These concepts are also common in other branches of 
condensed- matter physics [l6| . 

The entropy production constitutes the nonequilib- 
rium potential of the system which not far from equilib- 
rium allows us to infer regression laws relating currents 
and the conjugated thermodynamic forces. Thus in our 
case, from Eq. ^ we find 
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with L (p) being the matrix of Onsager coefficients which, 
as required for the second law, should be positive-definite. 
By combining here Eqs. (jl]) and ([6]) we obtain 
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where D(p) = fcsL/p is the matrix of diffusion coeffi- 
cients. Whence, Eq. ([7]) yields 
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Since the quanta are massless particles, there is no dif- 
fusion in physical or real space, J^; = 0. Therefore, from 
the set of equations ([5]) we obtain the diffusion current 
in momentum space 
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which defines t(p), the relaxation time simply related to 
the effective diffusion along the reaction coordinate x. 
Additionally, the size of an elementary cell in the phase 
space is given by h and, since according to the theory 
of Brownian motion the diffusion coefficient Z?, basically 
the mean square displacement (an area) per unit of time, 
is in direct relationship to the surface visited in the unit 
time in the phase space, wc can therefore interpret D as 
the ratio between the size of a characteristic cell, /i, and 
a characteristic time, t(p). 

By integrating Eq. ([9]) through from xi to 2:2, the 
position of both NPs, respectively, we find the net current 
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where pj{p,t) — p{p,x = j,t) is understood as related 
to the population of quasi-particles at Xj . Here, we have 
defined the net current as 
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where 



it) 



I pr{p)dr 



(12) 



stands for a hierarchy of relaxation times corresponding 
to a non-Deb ye p rocess proper of the dynamics of com- 
plex systems [l7j. 

3. Stationary state and thermal conductivity. — At 
equilibrium p2 = pi and J = 0. However, if we maintain 
the system in a stationary state in which each one of the 
NPs remains in local equilibrium with its respective bath 
at different temperatures so that pj pj^eq.(Tj) with 



eyi-p{-f3jEn) 
E„exp(-/3j£;„) 



(13) 



being the canonical distribution, where (3j = l/ksTj and 
En refers to the energy levels of the NP, then the current 
given through Eq. PH)) reaches a non-zero stationary 
value. 
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The fact that each nanoparticle remains in local equilib- 
rium with its proper bath corresponds with one among 
several possible metastable states of the composed sys- 
tem. 

We find the energy flux by multiplying Eq. (fT4|l 
through by En and canonically averaging, yielding 
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where we have defined the frequency uj — A/h, with A 
representing the gap between energy levels, 0(a;,T) — 
fiujN{u},T) is the average energy of an harmonic oscilla- 
tor and N{lo,T) = 1/ {expihuj /ksT) - 1) is the Planck 
distribution. The heat conductance follows by linearizing 
Eq. with respect to the temperature difference: 
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where Tq = (Ti -|-r2)/2 is the temperature corresponding 
to the stationary state of the system. It should be noted 
that these two spheres of the same diameter maintained 
at different temperatures constitute a model correspond- 
ing to more realistic systems such as for example a hot 
tip in contact with a flat substrate, as mentioned above. 
When d — 2'kc/uj, a phonon-like dispersion relation, Eq. 
(fT6)) gives us the heat conductance as a function of the 
distance between the NPs. Here, the effective relaxation 
time T*{uj) plays the role of an adjustable parameter 
which in general depends on the frequency according to 
the fact that for extremely close distances the system 
adopts a glassy behavior typical of complex systems. 
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FIG. 1: (Color online) Thermal conductance G12 vs dis- 
tance d reproducing the molecular dynamics data obtained 
by Domingues et al. The grey points represent the con- 
ductance when the particles with eflFective radius R = 0.72, 
1.10, and 1.79 nanometers are in contact. The lines show the 
analytical result obtained from Eq. (|16|l by adjusting r*{uj) 
to the simulation data. 



The intricate slow dynamics which complex systems 
exhibit is represented by a hierarchy of relaxation times 
often given by a Kohlrausch- Williams- Watts law [l^ 
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also known as stretched exponential and where /3 is a 
matching parameter and a and b are constants which de- 
pend on the size of the NPs. This is an empirical time- 
law, ubiquitous in disordered systems such as glasses 
which can be related to the distribution of residence times 
or the response function [13, [ill- Therefore, in our case 
we have tried a time law given by a stretched exponential 
to make Eq. (jl6p fit the available numerical simulation 
data [l[. We have found that an exponent /3 = 2 for 
which T*{t) as well as its Fourier transform t*{uj) are 
Gaussians, gives an excellent fit. This is supported by 
Fig. ([1]), where we have represented the heat conduc- 
tance as a function of the distance d between the NPs 
of different effective radii R — D/2. This figure shows 
a strong enhancement of the heat conductance when d 
decreases until around 2D nm due to multipolar inter- 
actions 0. When both NPs are in contact, a sharp fall 
occurs which can be interpreted as due to an intricate 
conglomerate of energy barriers inherent to the amor- 
phous character of these NPs generated by the strong 
interaction. The behavior shown in Fig. ^ agrees with 
what Eq. (flB]) predicts since in this equation, the two 
factors kBh/AT*{uj) and {{hukBTo)/ sinh{huj/2kBTo)y 
compete. For large distances the dominant factor is the 
latter factor whereas at contact the former factor be- 
comes the leading term. For large length scales, when r* 
is almost a constant which corresponds with exponential 
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relaxation the fluctuation-dissipation theorem is satis- 
fied and we find the system in the fluctuation-dissipation 
regime. Therefore, the conclusions of our theory are 
clearly well supported by Fig. ^ which indicates that 
the hypothesis of the rough energy landscape is indeed 
reasonable. 

As far as we know the decay in the thermal conductiv- 
ity at contact is an open question not yet resolved. Here, 
based on the fact that the crystalline structure of the 
NPs rapidly disappears when the intcrnanoparticle dis- 
tance is decreased the assumption that the resulting 
amorphous material self-organizes into an irregular distri- 
bution of energy barriers is the key to interpreting the ori- 
gin of this decay. Namely, the classical maximum thermal 
conductance applies to the fluctuation-dissipation regime 
before contact in which the system preserves its crys- 
talline structure and optical phonons are mainly excited. 
However, since it has been observed that both NPs be- 
come amorphous-like particles when they come closer to- 
gether due to their strong interaction, we have supposed 
that they behave as glassy materials. In glasses as well 
as in other disordered materials, it is a known fact that 
the distribution of modes presents anomalies which result 
from an excess of low- frequency modes [20] , the so-called 
boson peak. Hence, it is not possible to make a clear 
distinction between optical and acoustic modes. These 
low-frequency modes involve collective motions causing 
extremely slow structural relaxation. Therefore, this ac- 
counts for the high non-linear behavior of the thermal 
conductance between both NPs. 

4- Conclusion. — Up to now nanoscale radiative heat 
transfer has been studied assuming the validity of the 
fluctuation-dissipation theorem related to exponential or 
Debye relaxation. This approach shows itself to be suc- 
cessful to describe the energy exchange between two NPs 
before contact. At contact, however, since the interaction 
between both NPs does not admit a multipolar expansion 
this approach ceases to be valid. Hence, in view of the 
subsequent amorphous character of both NPs, we have 
assumed that this interaction gives rise to an irregular 
array of energy barriers and therefore, the relaxation of 
both NPs becomes a complex activated process. Thus, 
the transfer of energy is assumed to be due to the diffu- 
sion of quasi-particles through a rough energy landscape 
which gives rise to activated processes related to a rate 
current derivable in the framework of our thermodynamic 
theory. Our theory provides us the expression of the ther- 
mal conductance which depends on an adjustable param- 
eter representing a hierarchy of relaxation times, analo- 
gous to the slow relaxation behavior exhibited by similar 
complex systems. We have found that a time law similar 
to a stretched exponential with an exponent to the sec- 
ond power perfectly fits the available molecular dynamics 



simulation data, thus supporting our approach. 
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